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To elucidate the role in hadron structure of quark zero modes associated with instantons, the lowest eigen- 
functions of the Dirac operator have been calculated in quenched QCD, full QCD, and full QCD with cooling. 
Eigenmodes associated with a single instanton and an instanton-anti-instanton pair were studied to understand 
the role of lattice artifacts for Wilson Fermions. By truncating the spectral representation of the quark prop- 
agator, we show that the rho and pion contributions to the vector and pseudoscalar correlation functions are 
dominated by the zero-mode contributions. 



1. Introduction and Motivation 

This work is part of an ongoing effort to use 
the lattice to obtain insight into the structure 
of light hadrons by finding the paths that dom- 
inate the QCD path integral [Q. The picture 
that emerges is that of the instanton liquid de- 
veloped by Shuryak]^ and by Dyakanov and 
Petrov|^]. In this picture, light quarks propagate 
by hopping between zero modes, and the result- 
ing 't Hooft interaction]^ accounts for many phe- 
nomena, ranging from the channel dependence of 
vacuum hadron current correlation functions |^] 
to the so-called "spin crisis" in which helicity is 
transferred from valence quarks to gluons and qq 
excitations and the dd to uu ratio in the proton. 
Previous calculations show that removing non- 
instanton components of the gluon configurations 
by cooling reproduces the essential features of 
the instanton model and observables calculated 
with all gluons are similar to those obtained us- 
ing only instantons that survive cooling|^]. 

The chief limitation of these cooling studies 
is that although the loss of a single instanton 
by cooling can be avoided by the use of an im- 
proved action[^, the annihilation of instanton- 
anti-instanton pairs is unavoidable S. Hence, the 
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purpose of this work is to avoid the limitations 
of cooling by studying the low eigenfunctions of 
the Dirac operator directly in uncooled configu- 
rations. 

One essential piece of background information 
concerns vacuum matrix elements of meson cor- 
relation functions {ft\T J (x) J in the vec- 
tor and pseudoscalar channels 0. Lattice calcu- 
lations agree well with phenomenological disper- 
sion analyses]^, and the strength in the vicinity 
of 1 fm is strongly dominated by the rho and pion 
contributions. 

2. Eigenmodes of the Dirac operator 

Using the k-step Arnoldi method |jl^ we have 
calculated eigenmodes of the Wilson Dirac oper- 
ator 

/J 

+ (r + 7m)Wx-p,m^^- 

in several simple background configurations to 
understand their structure, as well as on 16^ lat- 
tices for quenched QCD at /3 = 5.85 and for un- 
quenched QCD at /3 = 5.5 and Kg^a = 0.16. 

As is well known, in the free discrete 
set of points starting at the lowest Matsubara 
mode (0, 0, 0, zbvr/Lf) approximates the contin- 
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Figure 1. Lowest 64 eigenvalues of the Dirac op- 
erator for an unquenched gluon configuration. 

uum spectrum ^[m + i|p|] in the physical regime 
and unphysical fermion modes are pushed to large 
(real) masses. In the presence of an instanton of 
size p at a; = 0, there are two lattice artifacts. Be- 
cause of the Wilson term, the eigenvalue acquires 
a shift of approximately i?eAo — (V'ol-C'IV'o) — 
Sku/p^. Also, in addition to a mode with zero 
imaginary part that approaches the continuum 
result 



tpo{x)s,a = Us 



V2 



TT (a;2 + p2)3/2 



one spin and color component of the lowest Mat- 
subara mode (0, 0, 0, ±7r/Lt) mixes with the zero 
mode with a coefficient that goes to zero in 
the large volume limit. The eigenvalues for an 
instanton-anti-instanton pair are approximately 



(A/ - Xa)' 



(Tia) 



2 V 4 

so that when the interaction Tja exceeds the dif- 
ference in the Wilson term shifts, complex conju- 
gate pairs of eigenvalues move slightly off the real 
axis. 

Fig. |l| shows the lowest 64 complex eigenvalues 
of the Dirac operator on an unquenched gluon 
configuration. Note that the eigenvalues of this 
uncooled configuration have just the structure 
expected from the pure one- and two-instanton 
studies above: isolated zeros along the real axis 
corresponding to individual instantons and com- 
plex conjugate pairs corresponding to instanton- 
anti-instanton pairs. To set the scale, note that 
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Figure 2. Contributions of low Dirac eigenmodes 
to the vector (upper graph) and pseudoscalar 
(lower graph) vacuum correlation functions. 

the lowest Matsubara mode would occur at 0.06 
on the imaginary axis, corresponding to 380 MeV, 
so all the conjugate pairs below this value are un- 
derstood to be the results of zero modes. 

3. Zero mode expansion 

The Wilson-Dirac operator has the prop- 
erty that D = 7515^75, which implies that 
("0^175 1 V'j) = unless Ai = A* and we may write 
the spectral representation of the propagator 

{x\il^i){'4^\i5\y) 



('0d75|^»)Aj 



where A,- = \* . A clear indication of the role 
of zero modes in light hadron observables is the 
degree to which truncation of the expansion to 
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the zero mode zone reproduces the result with 
the complete propagator. 

Fig. I shows the result of truncating the vac- 
uum correlation functions for the vector and 
pseudoscalar channels to include only low eigen- 
modes||ll[. On a 16^ lattice, the full propagator 
contains 786,432 modes. The top plot of Fig. || 
shows the result of including the lowest 16, 32, 
64, 96, and finally 128 modes for an unquenched 
configuration with a 63 MeV valence quark mass. 
Note that the first 64 modes reproduce most of 
the strength at 1 fm (10 lattice spacings) where 
the correlator is saturated by the rho resonance 
peak, and by the time we include the first 128 
modes, all the strength is accounted for. Simi- 
larly, the lower plot for a quenched configuration 
with a 23 MeV quark mass shows that the lowest 
128 modes also account for the analogous pion 
contribution to the pseudoscalar vacuum correla- 
tion function. Thus, without having to resort to 
cooling, by looking directly at the contribution of 
the lowest eigenfunctions, we have shown that the 
zero modes associated with instantons dominate 
the propagation of rho and pi mesons in the QCD 
vacuum. 

3.1. Localization 

Finally, it is interesting to ask whether the lat- 
tice zero mode eigenfunctions are localized on in- 
stantons. This was studied by plotting the quark 
density distribution for individual eigenmodes in 
the x-z plane for all values of y and t, and com- 
paring with analogous plots of the action den- 
sity. As expected, for a cooled configuration the 
eigenmodes correspond to linear combinations of 
localized zero modes at each of the instantons. 
What is truly remarkable, however, is that the 
eigenfunctions of the uncooled configurations also 
exhibit localized peaks at locations at which in- 
stantons are identified by cooling. Thus, in spite 
of the fluctuations several orders of magnitude 
larger than the instanton fields themselves, the 
light quarks essentially average out these fluctu- 
ations and produce localized peaks at the topo- 
logical excitations. When one analyzes a number 
of eigenfunctions, one flnds that all the instan- 
tons remaining after cooling correspond to local- 
ized quark fermion peaks in some eigenfunctions. 



However, some fermion peaks are present for the 
initial gluon configurations that do not coincide 
with instantons that survive cooling, correspond- 
ing to instanton-anti-instanton pairs that were 
annihilated during cooling. 
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